Explicit predictability and dispersion scaling 
exponents in fully developed turbulence 



Frangois G. Schmitt 

CNRS, UMR 8013 ELICO, Wimereux Marine Station, University of Lille 1 
28 av. Foch, 62930 Wimereux, France 



Abstract 

We apply a simple method to provide explicit expressions for different scaling expo- 
nents in intermittent fully developed turbulence, that before were only given through 
a Legendre transform. This includes predictability exponents for infinitesimal and 
non infinitesimal perturbations, Lagrangian velocity exponents, and dispersion ex- 
ponents. We obtain also new results concerning inverse statistics corresponding to 
exit-time moments. 
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1 Introduction 

One of the main feature of homogeneous and isotropic fully developed turbu- 
lence, corresponding to very large Reynolds numbers, is the fact that statistical 
properties at small scales are scaling and universal [1], and possess intermit- 
tent fluctuations [2]. This is classically characterized using structure functions 
scaling exponents, describing the probability density of velocity fluctuations 
at all scales belonging to the inertial range [3]. This concerns the Eulerian 
velocity and its statistics, but the characteristic features of fully developed 
turbulence is not only intermittency and multiple scales, it is also the loss of 
predictability, chaotic behaviour, dispersion and mixing properties, and a large 
number of degrees of freedom. In the last ten years, many of these properties 
have been revisited taking into account intermittency to consider how different 
scaling laws related to these issues are modified (for a review of several results 
in this framework see [4]). The scaling exponents obtained have often been ex- 
pressed as a Legendre transform of some function of velocity exponents, and 
numerical applications need a numerical Legendre transform involving the di- 
mension function of velocity fluctuations. Here we present a simple method to 
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simplify Legendre transform expressions and to replace them with paramet- 
ric functions that provide explicitely the scaling exponents when the Eulerian 
velocity scaling exponent structure functions are known. In the following the 
general framework is presented in section 2 and applied to different applica- 
tions from predictability to e-entropy studies in section 3. Let us note that 
the results presented here correspond to the ideal situation of homogeneous 
and isotropic turbulence; for realistic flows, this property cannot be expected 
to hold perfectly since the influence of boundary conditions may introduce 
statistical inhomogeneity. 



2 The multifractal framework and Legendre transform 

We recall first here the multifractal framework for velocity intermittency in 
turbulence, classically characterized using statistical moments of velocity in- 
crements. We recall the Legendre transform introduced to relate scaling mo- 
ment function to the codimension function of velocity singularities. The latter 
is subsequently systematically used to simplify the Legendre transform of some 
general expression of the codimension function. 

2. 1 Singularities, moment functions and Legendre transform 

In fully developed turbulence, intermittency is classically characterized by 

C{q), the scaling exponent of spatial structure functions. Denoting AVi = 
\V{x + i) — V{x)\ the increments of the velocity field at a spatial scale i, their 
fluctuations are characterized, in the inertial range, using the scale invariant 
moment function ({q) (see [3,5] for reviews): 

< AV^/ >~ (1) 

where q is the order of moment. Kolmogorov's initial proposal [1], for a non- 
intermittent constant dissipation, leads to C{q) — q/3. For intermittent turbu- 
lence, is a cumulant generating function, and is nonlinear and concave; 
only the third moment has no intermittency correction: ^(3) = 1 (see below). 
On the other hand, taking into account intermittency for the Eulerian velocity, 
one can write locally: 

AVe ~ (2) 

where /i is a singularity of the velocity fluctuations. There is a continuous range 
of possible h values, and they are characterized through their codimension 
function c{h) defined as [3,5]: 

p{AVe) ~ r^'^) (3) 
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where p{AVe) is the probabihty density of velocity fluctuations at scale i. 
Here the codimension is used for convenience, instead of the more classical 
dimension f{h) = d — c{h), where d is the dimension of the space {d = 1 
or d = 3 in most studies), and p(AV^) is the probability density of velocity 
fluctuations. Let us note that, while h can be both negative or positive (but 
it is most of the time positive), the codimension function c{h) is positive and 
decreasing. This continuous range of h values justifled the use of the term 
"multifractal" for such process [6]. 

The moments write: 

< AV,' > = j AV^piAV,) (4) 
~ j l''^^<^^p{h) ~ (5) 

using a saddle point argument [6] , this gives the classical Legendre transform, 
between Ce(?) the codimension function: 

C(?) = min{g/i + c(/i)} (6) 

n 

This can also be written in the following way, emphasizing the one-to-one 
relation between orders of moment q and singularities h: 

qh + c{h) = C{q) 

q = -c'{h) ^ ^ 

The last relation can be replaced by h = C{q)- Equations (7) are used below. 



2.2 Estimates of structure function scaling exponents up to moments of order 
8 



The scaling moment function C^{q) is proportional to the cumulant generating 
function of the variable log |AV^|; it is nonlinear and concave. There are two 
"flxed points" for this curve, corresponding to a piori exact values: C(0) = 
and C(3) = 1. The latter relation is a consequence of the Kolmogorov 4/5 
law [7] for the third order longitudinal structure function (without absolute 
values) : 

< AVl -Ul (8) 

The other moments do not have well-accepted theoretical values; the analytical 
curves depend on the stochastic model chosen for the statistics of velocity 
fluctuations. In agreement with Kolmogorov's initial hypothesis concerning 
the universality of velocity fluctuations in the inertial range, experimental 
estimates are rather stable for small and medium orders of moments and do 
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not depend on the Reynolds number: the exponents (e{<i) have been estimated 
experimentally for more than twenty years, for different types of flows (see 
[8,9,10,11,12,13]), and are now considered as rather stable and almost universal 
until moments of order 7 or 8 [10,11]. Let us note that the moment of order 
7 seems to be close to 2, for any reasons that did not receive up to now any 
theoretical founding. 

Since these estimates are rather stable up to moments of order 8, we chose in 
the following for simplicity for Ce{q) average values of the estimates published 
in Ref. [9,11,12,13]: see Table 1. These are estimates obtained from many 
experiments with Rx going from 35 to 5 000. The experimental average which 
is used here is quite close to all these experimental estimates. 



3 Applications to different predictability and dispersion turbu- 
lence functions 

3.1 Predictability of small perturbations in the intermediate dissipation range 
of turbulence 

Let us first consider the predictability property of small perturbations in tur- 
bulence. It is well-known that the viscous cutoff, corresponding to the scale 
for which the Reynolds number is of order 1, is different from the Kolmogorov 
scale rj when intermittency is taken into account: more intense singularities 
have a smaller local cutoff scale and weak fluctuations have larger cutoff scales. 
This range of scales is the intermediate dissipation range, where there is a mix- 
ture of inertial and dissipation effects. Using Eq. (2) and the Reynolds number 
Re = VL/p, this writes [14]: 

rjh = Li?e-i/(^+'^) (9) 

where L is a large scale and rjh is the cutoff scale associated to singularity 
h. This spatial scale is associated to a temporal cutoff scale through the 
velocity Vh = {rjh)^- Th = rjh/vh = {VhY^^ giving the singularity dependent 
cutoff time scale: 

Th = TRe'^ (10) 

where T is the large time scale associated to L. This temporal cutoff scale 
is the smallest scale of the system: below these scales viscosity dominates. 
In fully developed turbulence, the positive maximum Lyapunov exponent A, 
characterizing the exponential growth rate of an infinitesimal disturbance, is 
proportional to the inverse of the smallest characteristic time of the system. 
Since there is a range of time scales associated to the continuous range of h 
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values through Eq.(lO), the Lyapunov exponent may be given by [15,16]: 

A ~ y —dp{h) ~ J ReT^Re^dh (11) 

giving finally: 

A ~ Re" (12) 

with the exponent a given by a saddle-point argument since the Reynolds 
number is large [15,16]: 

(1-h- c(h) ] , , 

a = max <^ } 13 

h \ 1+h J ^ ^ 

This was solved numerically by Crisanti et al. [15] to provide a ~ 0.46 through 
a numerical Legendre transform using experimental estimates of the function 
c{h). This result is slightly smaller than the one provided by Ruelle [17] for 
nonintermittent turbulence: a — 1/2. 

Let us go further and show here how to obtain this exponent, considering, for 
more generahty, the moments of order q > oi the inverse times: 

< ^ >~ i?e°(«) (14) 

where, using as done above a saddle-point argument, the exponents a{q) are 
given by: 

(15) 

In general, the singularity ho for which this condition is met verifies G'{ho) — 
where 

= "'^-^j-^W (16) 

Estimating the derivative G'{ho) — and using Eq. (7) to introduce the mo- 
ment go associated to the singularity ho leads to: 

C(go) = 2q-qo (17) 

This provides a unique value of qo and hence of ho, for a given value of q. 
This equation, together with another use of the Legendre transform, gives the 
function a{q) — G{ho): 

a{q) = G{ho) = - - (C(go) - 

1 + ho 

= ?o - g 
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This gives finally the following parametric relation between a{q) and the mo- 
ment function C{q)'- 



The nonintermittent curve obtained for ({q) = q/3 is a{q) = q/2. The ex- 
plicit parametric relation (18) can be used to represent the curve a{q) for 
experimental estimates of C{q)' this is shown in Fig. 1 using the average (^{q) 
curve. One can see from Fig. 1 that up to moments of order 2.5 the linear 
non- intermittent curve is a very good approximation. The value of a{l) = a 
is smaller than 1/2 (close to 0.48) but the intermittent correction is very small. 
For q — 1 we have a — a{l) = qo — 1, where qo is given by C{Qo) = 2 — go- 
For q — 2 one obtains go = 3 since ({3) — 1, and hence a{2) = 1 is a fixed 
point, non affected by intermittency, recovering a result already given in Ref. 
[15]. One can also obtain the following estimates: «(3) ~ 1.56, q;(4) ~ 2.2, 
q;(4.5) ~ 2.5 (coming from ^(7) ~ 2), and a(5) ~ 2.85. 



3.2 Predictability of noninfinitesimal perturbations in turbulence 

The same approach can be used to express scaling exponents characterizing 
time scale statistics associated to noninfinitesimal perturbations. As proposed 
in Ref. [18] and developed in Ref. [4], and using the notations of the lat- 
ter reference, we consider a perturbation of size 6 of the velocity field, for 6 
belonging to the inertial range (hence the term noninfinitesimal). The time 
scale associated to the eddy of typical velocity S can also be considered as 
the decorrelation time associated to the perturbation. In this framework, the 
inverse time is 1/r ~ S/£ and since using the singularity h, one has i ~ 5^^'^, 
we obtain locally the relation between the time scale and the perturbation S: 



As before, this can be used to consider the moments of order g > of the 
inverse times, but here we will consider their scaling property as a power-law 
function of the velocity perturbation 6. Using the same type of argument as 
above, the moments of the inverse time are expressed as an integral over all 
singularities h an we consider all functions as power-law functions of 5 so that 
using a saddle-point argument, one has a scaling law: 



C(go) = 2g - go 
a(?) ^Qo-q 



(18) 




1 



(19) 



r 



< 1 >^ s-m 



(20) 
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where the exponents j3{q) arc given by the following relation (using the fact 
that the probability density writes p{h) ~ £^^'^^ ~ ^c(/i)//ij. 



where the negative sign has been introduced for convenience to have positive 
final numerical values for P{q) (see below). As done above for infinitesimal 
perturbations, this can be expressed explicitely introducing the singularity 
ho that maximises the expression into brackets. With the same procedure as 
above: involving a diffcrenciation and the introduction of the moment of order 
qo associated to ho, we obtain the final result as a parametric relation between 
/3(g) and C(g): 

C(go) = q ^22) 
Pil) =qo-q = C^~^\q) - q 

The exponent (3{q) is simply linked to the reciprocal of (, denoted C*'^^''- The 
nonintermittent curve obtained for ({q) = q/3 is (3{q) = 2q. Using Eq.(22) 
we can see that C(3) = 1 gives /3(1) = 2, which is a fixed point for /3(g), 
non affected by intermittency corrections [18]. Figure 2 shows the resulting 
function compared to the straight line of equation 2g. It is visible that /?(g) 
grows very fast. Indeed, the relation ^(7) ~ 2 gives the value /3(2) ~ 5. We 
can see that for g < 1, /3(g) < 2g and for g > 1, /3(g) > 2g. The result /3(1) = 2 
was checked using shell-model simulations [4]. 



3.3 Lagrangian velocity structure functions' in turbulence 



As an analogy with Kolmogorov's dimensional analysis in the Eulcrian frame- 
work. Landau [19] proposed in 1944 a 1/2 law for the temporal increments of 
the Lagrangian velocity AVr = \V{t + r) — V{t)\. This was later generahzed 
by Novikov, with a Lagrangian intermittency framework for the velocity [20] : 

< AV^^ >~ r^^('') (23) 

As for the Eulerian case, for a constant dissipation one obtains the "mean 
field" expression, neglecting intermittency: Cl(Q') = Q'/2- In this framework, 
the third order moment for the Eulerian velocity is analogous to the second 
order moment for the Lagrangian velocity: in case of intermittency Cl(?) is 
nonlinear and concave, and the non-intermittent function is valid only for 
g = 2:CL(2) = l. 

Recently some authors have proposed an hypothesis helping to relate Eulcrian 
and Lagrangian structure functions [21,22,23,24], following an earlier proposal 
by Borgas [25]. They consider the velocity advecting Lagrangian trajectories 
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as a superposition of different velocity contributions coming from different 
eddies having different characteristic times. After a time r the fastest eddies, 
of scale smaller than ^, are decorrelated so that at leading order, they assume 
[21,22,23]: 

AK ~ AC/^ (24) 
writing AC/^ ~ £/r and using AC/^ ~ l^, gives the time and space local cor- 
respondence r ~ Introducing this relation inside the integral of Eq.(7), 
one obtains for the moments of the Eulerian velocity increments: 

r qh+c(h) 

< MJ1 >- / T^=^dp{h) (25) 

Using a saddle point argument and using Eqs.(24) and (23), this gives the 
Lagrangian structure function as a Legendre transform [21]: 

The authors that explored this relation did not go further to express Cl(9) iii 
the general case [21,22,23]. 

In fact, as we showed elsewhere for a Lagrangian turbulence study [26], using 
an approach analogous to the ones above, one can obtain an explicit relation 
between Eulerian and Lagrangian structure functions' scaling exponents. We 
introduce the singularity ho that maximises the expression into brackets, and 
with the same procedure involving a differenciation and the introduction of 
the moment of order associated to Hq we obtain the final result (see also 
[26]): 

Cdq) = (Eiqo) ^27) 

q = qo- Ce(9o) 

This relation provides Cl(?) when the fonction Ce(?) is known, the second lign 
giving the link between q and go- We can check that for qo — 3, the second 
line gives q = 2 and the first one, Cl(2) = Ce(3) = 1- We can also consider the 
case qo — 7 and using Ce{7) — 2, the approximate result Cl(5) c::^ 1. 

Let us note that an analogous expression is provided in [25] for singulari- 
ties and (with our notations) codimension functions c(h) in his Appendix A 
(Equation (A 4)). Borgas used for this other types of arguments, including 
a Reynolds number scaling expression for Eulerian and Lagragian statistics; 
however, he did not go to the moment framework to provide an equivalent 
expression for structure functions scaling exponents. Let us also note that 
some comparisons are performed in [21] between their shell model Lagrangian 
values and Eulerian shell model estimates, transformed through a numerical 
Legendre transform using the analytical expression of a fit of c{h). The same 
type of comparison is performed in [23] between their DNS estimates and a 
numerical Legendre transform of Eulerian values. The approach proposed here 
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is more direct since we test an explicit parametric relation corresponding to 
different hypothesis. 

The Eulcr- Lagrange relation given here is non-linear, so that one cannot ex- 
pect to obtain the same type of statistical models in Eulerian and Lagrangian 
frameworks. For example, as explained by Chevillard [24] using Borgas' rela- 
tion, in this framework, the statistics cannot be lognormal for both Eulerian 
and Lagrangian velocity fields. Indeed, for C^;(?) = aq — bq^ with a — {2 + ii)/6 
and b = /x/18, one obtains Ciiq) = (a - 1 + ^J(l^^aj^^^4qb - 2bq)/2b. For 
more details see our paper [26] where different proposals for the Lagrangian 
structure functions scaling exponents are tested. 

In Figure 3, the Lagrangian prediciton is compared to Lagrangian estimates 
published in recent papers [21,27,22]. These papers report different types of 
such estimates: Ref. [27] gives the values of experimental estimates for four 
flows going from Rx of 510 to 1 000, and Direct Numerical Simulations (DNS) 
estimates for Rx—75 and 140. Ref. [21] reports shell model estimates and Ref. 
[23] DNS corresponding to -Ra=280. In order not to have a too heavy figure, we 
plotted in Figure 3 four different types of Lagrangian values: an average value 
of the 4 experiments shown in Ref. [27], the DNS simulation corresponding to 
Rx—IAO of the same publication, the shell model estimates of Ref. [21] and 
the DNS estimate of Ref. [23]. The average Eulerian curve is transformed and 
compared to these four Lagragian estimates. Several comments can be pro- 
posed from Figure 3. First, let us note that these experimental and numerical 
estimations for Eulerian and Lagrangian scaling exponents come from different 
types of flows with different Reynolds numbers. For Eulerian estimates, the 
exponents have been carefully compared and are now considered as rather sta- 
bles up to moments of order 8, but for Lagrangian estimates, this work is still 
to do; measurements are only available for a few years, and the scaling ranges 
are still small, so that the range of variability for these exponents is expected 
to narrow in the future. Indeed, as for the Eulerian case, one can expect that 
for high Reynolds numbers these exponents are rather stable and universal, 
up to a given order of moment that may be smaller than the Eulerian case, 
for which this critical order of moments is about 7-8. 

With these comments in mind, we can provide some preliminary conclusions 
from the observation of Figure 3: (i) up to moment of order 3, the prediction 
are close to all Lagrangian data; (ii) for larger moments, there is too much 
scatter in the estimates to conclude about the validity of the transformation 
given by Eq.(27). A possible explanation of such scatter is the fact that scaling 
can be strongly affected by trapping in vortex filaments [23]. In this case, 
after a filtering of high frequencies, the data may be closer to the multifractal 
prediction. 
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3-4 Relative turbulent dispersion 



We next consider the influence of intermittency on relative dispersion: this 
deals with the statistics of particle pair dispersion, advected by an intermit- 
tent Eulerian velocity field. This topic is close to the Lagrangian case, let us 
consider R{t), the distance between a pair of particles advected by an inter- 
mittent turbulence. Following Ref. [20,28,29], one can introduce AVt-{R), the 
velocity difference along the Lagrangian trajectories. This can be written: 

uR 

^^AV, (28) 

where AVr is the increment of the Lagrangian velocity is evaluated for a 
distance R associated to a time scale r = R/AVr. Then using the same 
hypothesis as for the Lagrangian case (Eqs. (23) and (24)), we obtain the 
time evolution of the moments of the pair distance: 

< R" >-< {rAVrY >- r^(^) (29) 

where H{q) is the scaling exponent characterizing the intermittent pair disper- 
sion. Using Eq.(26) the exponents H{q) can be given as a Legendre transform 
of the Eulerian velocity, as provided by Boffetta et al. [28,29]: 

In fact Eq.(29) shows also directly that the dispersion exponent is simply 
related to the Lagrangian velocity exponent: 

H{q)^q + CL{q) (31) 
Equations (27) and (31) then give: 

(32) 

g = 90 - CE{qo) 

This can also be written H{q — C{q)) = q. This relation could also have been 
obtained from Eq.(30) using the same approach as above. With go = 3 we 
obtain q = 2 and H{2) = 3, corresponding to Richardson's law < R^ >~ 
[30], which is thus a fixed point for which there is no intermittency correc- 
tion, as already noticed [20,28,29]. In Fig. 4, we represent the transform (32) 
compared to the non-intermittent line of equation 3g/2, compared also to the 
shell model estimates published in Ref. [29]. This shows that the nonlinearity 
of H{q) does not seem very strong. Furthermore, H is concave, and the values 
give by Eq. (32) are quite close to the ones reported in Ref. [29] . 
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3.5 Exit-time moments and e- entropy 



We finally consider exit-time moments [31], also called inverse structure func- 
tions [32]. We consider a time series of a turbulent quantity X(t) having in- 
termittent fluctuations characterized by singularities h such as AX^ ~ 
with codimension c{h): p{h) ~ r'^W. The moments arc assumed to scale as 
< AX^ >~ r''*^'^^ with Eq.(9). One then studies the dynamics of such signal 
by considering exit times (also called distance structure functions or inverse 
turbulence functions): considering a threshold value 5 = AX^-, let us denote 
t{6) the flrst time needed to obtain an increment S = AXr- A new time series 
(tj) can be obtained this way, whose statistics will scale as the threshold value 
6: 

< >~ (33) 

where x{q) can be obtained, as before, as a Legendre transform involving c{h). 
Writing r ~ 5^^^ and inserting this into the integral of the moment estimation 
gives [32,31]: 

X(,) = miu(i±i<^) (34) 

Here as before, we estimate this function explicitely. A derivative of G{h) ~ 
{q + c{h))/h leads to x(g) = G(/to) with C(9o) = —(!■ We then obtain 

C(?o) = -q 
The result then writes simply: 

x(-C(go)) = -go (36) 

which is exactly a result given in Ref. [33] based on an exact result for a special 
case (a multifractal Cantor set), and experimentally verified using shell model 
turbulence simulations, but with no formal proof of the formulae in the general 
case. A similar expression obtained in the very different context of Laplacian 
random walk is given in Ref. [34]. 

We must note here, as already noticed in [33], that Eq.(36) involves for either 
the original series or for the return time series, negative moments. Let us con- 
sider flrst positive moments of return times: g > in Eq.(33). We need here 
the hypothesis of the existence of negative moments of the velocity structure 
functions, but it may not be the case, since this needs in practice extremely 
precise measurements. It is argued in [35] that negative moments for g < — 1 
do not exist since the pdf of velocity increments does not vanish at 0. The 
derivative writes G"(^o) = "(^ + C(Q'o))/^o- -P'o^ g > 0, if negative moments 
go such that C(5o) = — g < are not accessible in the experimental time se- 
ries, then the derivative stays strictly negative, G{h) does not reach a relative 
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minimum, and the minimal value selected by Eq.(34) will be obtained for the 
maximal value of h, which is a "minimal" singularity of the velocity fluctua- 
tions, and hence denoted here /imin: x{q) = G^(^min)- Since this value of h is 
the same for all g's, the result is a hnear behaviour: 

X{q) = (37) 

'''min 

The singularity h„im is associated to small moments and the smallest fluc- 
tuations detectable in experimental time series; it is not clear whether this 
maximum singularity reaches a flxed value for Eulerian turbulence, and what 
is this value. This maximal singularity also depends on the precision of the 
experimental measurements which are analyzed. To check this, we consider nu- 
merical and experimental analysis of such inverse statistics for g > 0. We plot 
in Fig. 5 the data for x(g) with g > reported in Refs. [32,36]. The latter are 
experimental estimates, while the former are obtained from shell model calcu- 
lations. Both are rather linear (the published experimental estimates contain 
only the flrst three moments) , and seem to conflrm a linear behaviour such as 
Eq.(37). Recently inverse statistics have been tested using wind-tunnel exper- 
imental data, conflrming a linear behaviour for positive (relative) exponents 
[35]. 

Let us now focus on negative moments of the time statistics: g > in Eq.(33). 
This corresponds to positive moments of inverse times. We have C(?o) — —Q 
and x{q) = G^ho) = —go- The result will be simply written introducing p = — g 
and go(p) = -x(g): 

< Q-J >~ (^-«°^) (38) 

with: 

CM=P (39) 

The exponent obtained with the inverse statistics is in fact the inverse (in the 
sense of reciprocal) of the structure function: 

qoip) = C^-'\p) (40) 

For Eulerian velocity (with Taylor's hypothesis to interprete time statistics 
as being close to spatial statistics), a time series measured at a flxed location 
will be characterized by ^ = (^e- Thus q = 1 selects go = 3, non affected by 
intermittency corrections. The non-intermittent curve is go = 3g, and taking 
into account intermittency we have also the approximate value go (2) ~ 7. The 
corresponding curves for Eulerian turbulence and for passive scalar turbulence 
(data from Rcf. [37] corresponding to an average of several experimental esti- 
mates) are shown in Fig. 6: they are increasing fast, especially for the passive 
scalar case. Let us note also that the experimental results published in [35] do 
not conflrm Eq.(40), except ior p — 1. 
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Exit-time statistics have also been used in Ref. [38,39,4] to characterize the 
e-entropy of turbulent signals. The concept of e-entropy has been introduced 
to quantify the predictability and degree of complexity of continuously vary- 
ing signals such as turbulent signals. The continuous signal is transformed 
into a symbolic time series using a grid in the phase space, of resolution e. 
The classical Shannon entropy h of the resulting symbolic time series is then 
studied, as a function of the grid size e (see [40,41]). Abel et al. [38,39] have 
proposed to study the e-entropy of turbulent signals h{e) using the exit-times 
t{e), the time for the signal to undergo a fluctuation of size e. They obtained 
the following expresion 

h{e) = h^{e) <-> (41) 
r 

where h^{e) is the exit-time e-entropy, being bounded and having only a log- 
arithmic variation with e. The leading e behaviour of h{e) is then given by 
< ^ > , where the average is done here by considering a succession of exit time 
events, and is thus related to Eq. (38). This gives using Eq.(38) and e — S: 

h(e) ~ e-^ (42) 

with P = qo{l) as given by Ref. [38,39]. But we can write also using Eq.(40): 

P^C^-'\1) (43) 

Hence for a turbulent time series recorded at a fixed location, but using Tay- 
lor's hypothesis, corresponding to Eulerian turbulence, P — 3 [38,39,4]. On the 
other hand, for Eulerian temperature turbulence (still with Taylor's hypothe- 
sis), one will have a larger value since Ce(3) — 0.8; using the values published 
in Ref. [37], we have C6'(4) — 1 and hence we obtain for passive scalar Eulerian 
turbulence I3g ~ 4. For times series of Lagrangian velocity and passive scalar 
turbulence, both characterized by ({2) = 1, we have (3 — 2. This recovers the 
result obtained in Ref. [40] , which have been obtained using a dimensional ar- 
gument, thus corresponding implicitely to Lagrangian velocity. Let us also note 
that Eq.(43) recovers the results of Kolmogorov [42]: for Brownian processes 
with scaling power spectra of the form E{uj) ~ uj~'^^^'^^\ he gave jS = 1/H. 
Equation (43) gives the same result, since for such process, ({q) = qH, hence 
(5 = C(-i)(l) = l/H. 



Conclusion 

We have introduced a simple procedure helping to express explicitely scaling 
exponents expressed as Legendre transforms in multifractal turbulence. We 
have applied this idea to several problems involving intermittent corrections 
to predictability and dispersion studies. Even if the basic idea is simple, this 
produced several interesting new results: several functions are obtained as 
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parametric transform of the Eulerian velocity structure function. We have 
confirmed and generahzed some previously known results and provided some 
predictions. 

More precisely, in this paper wc have expressed the scaling moment function 
of infinitesimal perturbations, that characterize the inverse times statistics in 
the intermediate dissipation range, as function of the Reynolds number. For 
finite size perturbations, we have also provided the scaling moment function 
P{q) of inverse times, as function of the size 5 of perturbations, showing that 
this exponent is very simply related to the reciprocal of the velocity struc- 
ture functions exponent ({q). We veryfied the known result j3{l) = 2 and 
obtain the approximate value /3(2) ~ 5. We have also considered Lagrangian 
velocity structure functions Cl(?) ^^'^ showed how to relate them explicitely 
to the Eulerian curve Ce{(i)- We compared the predicted Lagrangian curve 
to recent experimental values. We noticed the approximate value Cl(5) — 1. 
We also considered the scaling exponent H{q) for pair dispersion which is 
very simply related to the Lagrangian scaling exponents. The obtained curve 
was compared to recent numerical estimates. We finally considered exit-time 
moments, and confirmed in the general case a result obtained recently on ex- 
perimental grounds and analytically for a special case. We showed that, for 
exit-time positive moments, the scale invariant function may often be 
linear, due to finite precision of experimental measurements. This is in agree- 
ment with several published experimental results. We then proposed to focus 
on inverse statistics considering negative moments, hence positive moments of 
inverse times. In this framework we showed that the resulting exponent qo{q) 
is the reciprocal of C{q)- Using some recent results linking exist time statis- 
tics to e-entropy expression, we obtain a new general result for the e-entropy 
exponent of multi-affine signals (3 = ^'^"^•'(1) confirming (3 = 3 for Eulerian 
turbulence, and giving as prediction /3 = 2 for Lagrangian turbulence (velocity 
and passive scalar), /3 ~ 2 for financial time series, and the approximate value 
/3 ~ 4 for passive scalar Eulerian turbulence. 



Acknowledgements 

Useful suggestions by the referees are acknowledged. 



References 

[1] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 299 (1941). 

[2] G. K. Batchelor and A. A. Townsend, Proc. Roy. Soc. A 199, 238 (1949). 



14 



[3] U. Frisch, Turbulence; the legacy of A. N. Kolmogorov (Cambridge University 
Press, Cambridge, 1995) 296 pp. 

[4] G. Boffctta and M. Cencini and M. Falcioni and A. Vulpiani, Phys. Rep. 356, 
367 (2002). 

[5] D. Schertzer and S. Lovejoy and F. Schmitt and Y. Chigirinskaya and D. 
Marsan, Fractals 5, 427 (1997). 

[6] G. Parisi and U. Frisch, in Turbulence and Predictability in Geophysical Fluid 
Dynamics and Climate Dynamics, M. Ghil and R. Benzi and G. Parisi eds.. 
North Holland (Amsterdam, 1985) 84. 

[7] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 32, 19 (1941). 

[8] F. Anselmet and Y. Gagne and E. J. Hopfinger and R. A. Antonia, J. Fluid 
Mech. 140, 63 (1984). 

[9] R. Benzi et al., Phys. Rev. E 48, R29 (1993). 

[10] F. Schmitt and D. Schertzer and S. Lovejoy and Y. Brunet, Nonlinear Proc. 
Geophys. 1, 95 (1994). 

[11] A. Arncodo et al., Europhys. Lett. 34, 411 (1996). 

[12] R. A. Antonia and B. R. Pearson, Europhys. Lett. 40, 123 (1997). 

[13] W. van de Water and J. A. Herweijer, J. Fluid Mech. 387, 3 (1999). 

[14] G. Paladin and A. Vulpiani, Phys. Rep. 156, 147 (1987). 

[15] A. Crisanti and M. H. Jensen and A. Vulpiani and G. Paladin, Phys. Rev. Lett. 
70, 166 (1993). 

[16] A. Crisanti and M. H. Jensen and G. Paladin and A. Vulpiani, J. Phys. A: 
Math. Gen. 26, 6943 (1993). 

[17] D. Ruelle, Phys. Lett. A 72, 81 (1979). 

[18] E. Aurell and G. Boffetta and A. Crisanti and G. Paladin and A. Vulpiani, 
Phys. Rev. Lett. 77, 1262 (1996). 

[19] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (MIR, First Russian Edition, 
1944). 

[20] E. A. Novikov, Phys. Fluids A 1, 326 (1989). 

[21] G. Boffetta and F. De Lillo and S. Musacchio, Phys. Rev. E 66, 066307 (2002). 

[22] L. Biferale and G. Boffetta and A. Celani and B. J. Devenish and A. Lanotte 
and F. Toschi, Phys. Rev. Lett. 93, 064502 (2004). 

[23] L. Biferale and G. Boffctta and A. Celani and A. Lanotte and F. Toschi, Phys. 
Fluids 17, 021701 (2005). 



15 



[24] L. Chevillard and S. G. Roux and E. Levcque and N. Mordant and J.-F. Pinton 
and A. Arneodo, Phys. Rev. Lett. 91, 214502 (2003). 

[25] M. S. Borgas, Phil. Trans. R. Soc. Lond. A 342, 379 (1993). 

[26] F. G. Schmitt, Phys. Rev. E (submitted). 

[27] N. Mordant and E. Leveque and J.-F. Pinton, New Journal of Physics 6, 116 
(2004). 

[28] G. Boffetta and A. Celani and A. Crisanti and A. Vulpiani, Europliys. Lett. 46, 
177 (1999). 

[29] G. Boffetta and A. Celani and A. Crisanti and A. Vulpiani, Phys. Rev. E 60, 
6734 (1999). 

[30] L. F. Richardson, Proc. Roy. Soc. London A 110, 709 (1926). 

[31] L. Bifcralc and M. Cencini and D. Vergni and A. Vulpiani, Phys. Rev. E 60, 

6295 (1999). 

[32] M. H. Jensen, Phys. Rev. Lett. 83, 76 (1999). 

[33] S. Roux and M. H. Jensen, Phys. Rev. E 69, 016309 (2004). 

[34] M. B. Hastings, Phys. Rev. Lett. 88, 055506 (2002). 

[35] B. R. Pearson and W. van de Water Phys. Rev. E 71, 036303 (2005). 

[36] S. Beaulac and L. Mydlarski, Phys. Fluids 16, 2126 (2004). 

[37] F. G. Schmitt, Eur. Phys. J. B (submitted). 

[38] M. Abel et al., Physica D 147, 12 (2000). 

[39] M. Abel et al., Phys. Rev. Lett. 84, 6002 (2000). 

[40] X. J. Wang and P. Gaspard, Phys. Rev. A 46, 3000 (1992). 

[41] P. Gaspard and X. J. Wang, Phys. Rep. 235, 291 (1993). 

[42] A. N. Kolmogorov, IRE Trans. Inf. Theory 1, 102 (1956). 



16 



Table 1 

Some recent experimental estimations for CeIq) and the average value used here. 
Ref. [9]: Re from 6,000 to 47,000; Ref. [11]: various experiments, Rx going from 
35 to 5,000; Ref. [12]: turbulent wake, R\=230; Ref. [13]: various experiments, R\ 



;oin] 


g from 340 to 800. 






q 


[9] 


[11] 


[12] 


[13] 


Average ^ 


1 








0.368 


0.37 


2 


0.71 


0.70 


0.69 


0.694 


0.70 


3 


1 


1 


1 


1 


1 


4 


1.28 


1.25 


1.29 


1.282 


1.28 


5 


1.53 


1.50 


1.55 


1.541 


1.53 


6 


1.78 


1.75 


1.79 


1.782 


1.78 


7 


2.01 


1.97 


2.01 


2.007 


2.00 


8 


2.22 


2.10 


2.22 


2.217 


2.19 
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1 2 3 4 5 

q 



Fig. 1. The curve a(g) (open dots: experimental values; continuous line: lognormal 
fit) compared to the non-intermittent curve of equation q/2 (dotted line). It is visible 
that a(l) = a < 1/2, but the intermittency correction is very small. The value 
a(2) = 1 is a fixed point not affected by intermittency. One obtains numerically 
a(3) ~ 1.56 and a(4) ~ 2.2. 
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0.5 1 1.5 2 2.5 



q 

Fig. 2. The curve P{q) (open dots: experimental values; continuous line: lognormal 
fit) compared to the non- intermittent curve of equation 2q (dotted line). This curve 
is strongly nonlinear and increasing very fast. The value /3(1) = 2 is a fixed point 
and we may note also that /3(2) ~ 5. 
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Fig. 3. The curve Cl(9) (open dots: experimental values; continuous line: lognormal 
fit) compared to the non-intermittent curve of equation g/2 (dotted line), and to 
four Lagrangian recent estimates: average of 4 experiments published in [27]; DNS 
from [27], with ^2^=140; shell model results from [21]; and DNS from [23]. The 
agreement is good for low orders of moment; for larger moments experimental and 
numerical estimates of scaling exponents have a larger scatter, and it is difficult to 
conclude concerning the validity of Eq. (27). 
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1 2 3 4 5 6 

q 

Fig. 4. The curve H{q) (open dots: experimental values; continuous line: lognormal 
fit) compared to the non- intermittent curve of equation 3q/2 (dotted line). This 
is also compared to the shell- model estimates published in [29]. The agreement is 
quite good. The curve H is only slightly nonlinear and concave. 
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Fig. 5. Representation of data for x{q) for g > reported in Refs. [32,36]. These 
correspond to numerical estimates from a shell model [32] and experimental esti- 
mates with Rx = 582 for velocity (longitudinal structure functions) and temperature 
turbulence [36]. These values are compatible with a linear behaviour (linear fits are 
displayed as dotted lines), as predicted by Eq.(37), although experimental estimates 
are provided here for only three values. 
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Fig. 6. The curve qoip) = C^~^\p) representing the scahng exponent function for 
the moments of order p > of the inverse exit times, for the velocity (open dots: ex- 
perimental values; continuous line: lognormal fit) compared to the non-intermittent 
curve of equation 3q (dotted line); also shown: the same curve for passive scalar tur- 
bulence (continuous line and closed dots, data from Ref. [37]). The passive scalar 
curve is very close to the velocity one for weak moments; it increases very quickly. 
We also represent the e-entropy exponent /? = to emphasize the clear 

different values for velocity and passive scalar Eulerian turbulence 
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